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Utilizing various gauges of the radial coordinate we give a description of static spherically sym- 
metric space-times with point singularity at the center and vacuum outside the singularity. We show 
that in general relativity (GR) there exist a two-parameters family of such solutions to the Einstein 
equations which are physically distinguishable but only some of them describe the gravitational 
field of a single massive point particle with nonzero bare mass Mo. In particular, we show that 
the widespread Hilbert's form of Schwarzschild solution, which depends only on the Keplerian mass 
M < Mo, does not solve the Einstein equations with a massive point particle's stress-energy tensor 
as a source. Novel normal coordinates for the field and a new physical class of gauges are proposed, 
in this way achieving a correct description of a point mass source in GR. We also introduce a gravi- 
tational mass defect of a point particle and determine the dependence of the solutions on this mass 
defect. The result can be described as a change of the Newton potential Lp N — —G N M/r to a modi- 
fied one: <p G = —G N M / (r + G N M/c 2 In ^j-) and a corresponding modification of the four-interval. 
In addition we give invariant characteristics of the physically and geometrically different classes of 
spherically symmetric static space-times created by one point mass. These space-times are analytic 
manifolds with a definite singularity at the place of the matter particle. 
PACS number(s): 04.20.Cv, 04.20. Jb, 04.20.Dw 



I. INTRODUCTION 

The Einstein's equations: 

G$ = kT? (1) 

determine the solution of a given physical problem up 
to four arbitrary functions, i.e., up to a choice of coor- 
dinates. This reflects the well known fact that GR is a 
gauge theory. 

According to the standard textbooks the fixing of 
the gauge in GR in a holonomic frame is represented by 
a proper choice of the quantities 




which emerge when one expresses the 4D d'Alembert op- 
erator in the form 5^V M V„ = (d M d„ - f^A). 
We shall call the change of the gauge fixing expressions 
, without any preliminary conditions on the analytical 
behavior of the used functions, a gauge transformations 
in a broad sense. This way we essentially expand the class 
of admissible gauge transformations in GR. Of course this 
will alter some of the well known mathematical results 
in the commonly used mathematical scheme of GR. We 
think that a careful analysis of such wider framework 
for the gauge transformations in GR can help us clarify 
some long standing physical problems. After all, in the 



*E-mail: fiziev@phys.uni-sofia.bg 



physical applications the mathematical constructions are 
to reflect in an adequate way the properties of the real 
physical objects. 

It is well known that in the gauge theories we may 
have different solutions with the same symmetry in the 
base space. Such solutions belong to different gauge sec- 
tors, owning different geometrical, topological and phys- 
ical properties. For a given solution there exists a class 
of regular gauges, which alter the form of the solution 
without changing these essential properties. In contrast, 
by performing a singular gauge transformation one can 
change both the geometrical, topological and physical 
properties of the solution, making a transition to another 
gauge sector. 

In geometrical sense the different solutions in GR de- 
fine different 4D pseudo-Ricmannian space-time mani- 
folds A4^ 1 ' 3 ^{g tll/ (x)}. A subtle point in GR formalism 
is that transitions from a given physical solution to an 
essentially different one, can be represented as a "change 
of coordinates" . This is possible because from gauge the- 
ory point of view the choice of space-time coordinates in 
GR means simultaneously two different things: choosing 
a gauge sector in which the solution lives and at the same 
time fixing the (regular) gauge in this sector. While the 
latter solves an inessential local gauge fixing problem, the 
first fixes the essential global properties of the solution. 

One can understand better this peculiarity of GR in 
the framework of its modern differential geometrical de- 
scription as a gauge theory, using the principal frame 
bundle F {M^ 1 '^ {g pLy {x)}) . (See, for example, the sec- 
ond edition of the monograph @], Q, and the refer- 
ences therein.) The change of coordinates on the base 
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space M^ 1 ' 3 ^ {g fJll y(x)} induces automatically a nontrivial 
change of the frames on the fibre of frames. Singular 
coordinate transformations may produce a change of the 
gauge sector of the solution, because they may change 
the topology of the frame bundle, adding new singu- 
lar points and(or) singular sub-manifolds, or removing 
some of the existing ones. For example, such transforma- 
tions can alter the fundamental group of the base space 
Ai^ 1 ' 3 ^ , {g fiu (x)}, the holonomy group of the affine con- 
nection on F (^M^ 1 ' 3 ^ {g fJ-u (x)}) , etc. Thus one sees that 
coordinate changes in broad sense are more than a pure 
alteration of the labels of space-time points. 

If one works in the framework of the theory of smooth 
real manifolds, ignoring the analytical properties of the 
solutions in the complex domain, one is generally allowed 
to change the gauge without an alteration of the phys- 
ical problem in the real domain, i.e. without change of 
the boundary conditions, as well as without introduc- 
tion of new singular points, or change of the character 
of the existing ones. Such special type of regular gauge 
transformations in GR describe the diffcomorphisms of 
the real manifold M.^ 1,3 ^ 1 {g^vix)}. This manifold is al- 
ready fixed by the initial choice of the gauge. Hence, the 
real manifold M^ 1 ' 3 ^ {g liy {x)} is actually described by a 
class of equivalent gauges, which correspond to all diffeo- 
morphisms of this manifold and are related with regular 
gauge transformations. 

The transitions between some specific real manifolds 
■W' 3 '{"Smn( r )}i which are not diffeomorphic, can be 
produced by the use of proper singular gauge transforma- 
tions. If we start from an everywhere smooth manifold, 
after a singular transformation we will have a new man- 
ifold with some singularities, which are to describe real 
physical phenomena in the problem at hand. The inclu- 
sion of such singular transformations in our consideration 
yields the necessity to talk about gauge transformations 
in a broad sense. They are excluded from present-day 
standard considerations by the commonly used assump- 
tion, that in GR one has to allow only diffeomorphic map- 
pings. 

Similar singular transformations are well known in 
gauge theories of other fundamental physical interac- 
tions: electromagnetic interactions, electroweak interac- 
tions, chromodynamics. For example, in the gauge theo- 
ries singular gauge transformations describe transitions 
between solutions in topologically different gauge sec- 
tors. Singular gauge transformations are used in the the- 
ory of Dirac monopole, vortex solutions, t'Hooft-Poliakov 
monopoles, Yang-Mills instantons, etc. See, for example, 
0, Q and the references therein. 

The simplest example is the singular gauge transfor- 
mation of the 3D vector potential: A(r) — > A(r) + 
Vtp(r) in electrodynamics, defined in Cartesian coor- 
dinates {x, y, z} by the singular gauge function ip — 
a arctan(y/x) (a — const). Suppose that before the 
transformation we have had a 3D space g mn (r) = 

} = TZ^{S mn }. Then this singular gauge transfor- 
mation removes the whole axes OZ out of the Euclidean 



3D space 71^ {8^}, changing the topology of this part of 
the base space. As a result the quantity <f c drA(r), which 
is gauge invariant under regular gauge transformations, 
now changes its value to § c drA(r) + 2irNa, where N is 
the winding number of the cycle C around the axes OZ. 
Under such singular gauge transformation the solution of 
some initial physical problem will be transformed onto a 
solution of a completely different problem. 

At present the role of singular gauge transformations 
in the above physical theories is well understood. 

In contrast, we still don't have systematic study of the 
classes of physically, geometrically and topologically dif- 
ferent solutions in GR, created by singular gauge trans- 
formations, even in the simple case of static spherically 
symmetric space-times with only one point singularity, 
although the first solution of this type was discovered 
first by Schwarzschild nearly 90 years ago Q . Moreover, 
in GR at present there is no clear understanding both of 
the above gauge problem and of its physical significance. 

Here we present some initial steps toward the clarifica- 
tion of the role of different GR gauges in broad sense for 
spherically symmetric static space-times with point sin- 
gularity at the center of symmetry and vacuum outside 
this singularity. 



II. STATIC SPHERICALLY SYMMETRIC 
SPACE-TIMES WITH POINT SOURCE OF 
GRAVITY 

The static point particle with bare mechanical rest 
mass Mo can be treated as a 3D entity. Its proper frame 
of reference is most suitable for description of the static 
space-time with this single particle in it. We prefer to 
present the problem of a point source of gravity in GR 
as an ID mathematical problem, considering the depen- 
dence of the corresponding functions on the only essential 
variable - the radial variable r. This can be achieved in 
the following way. 

The spherical symmetry of the 3D space reflects ad- 
equately the point character of the source of grav- 
ity. A real spherically symmetric 3D Riemannian space 
M^i-gmnir)} C M ( - 1 '^{g liu (x)} can be described 
using standard spherical coordinates r € [0, oo), 9 € 
[O,tt),0 € [0,2tt). Then 

r = {a; 1 , x 2 , x 3 } = {r sin 9 cos (f>,r sin 9 sin (j),r cos 9}. 

The physical and geometrical meaning of the radial 
coordinate r is not defined by the spherical symmetry of 
the problem and is unknown a priori [|| . The only clear 
thing is that its value r — corresponds to the center of 
the symmetry, where one must place the physical source 
of the gravitational field. In the present article we assume 
that there do not exist other sources of gravitational field 
outside the center of symmetry. 

There exists unambiguous choice of a global time t due 
to the requirement to use a static metric. In proper units 
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it yields the familiar form of the space-time interval |9(: 

ds 2 =gtt(r)dt 2 +g rr (r)dr 2 -p(r) 2 (d6 2 + sin 2 d<t> 2 ) (3) 

with unknown functions gtt{f) > 0, g r r{f) < 0, p(r). 

In contrast to the variable r, the quantity p has a clear 
geometrical and physical meaning: It is well known that 
p defines the area A p = Anp 2 of a centered at r = 
sphere with " a radius" p and the length of a big circle on 
it l p = 2np. One can refer to this quantity as an "area 
radius", or as an optical "luminosity distance", because 
the luminosity of distant physical objects is reciprocal to 

One has to stress that in the proper frame of refer- 
ence of the point particle this quantity does not mea- 
sure the real geometrical distances in the corresponding 
curved space-time. In contrast, if p fixed is some arbi- 
trary fixed value of the luminosity distance, the expres- 

1/2 

sion (p2 — pi) I (1 — 2M/pf ixe d) 1 measures the 3D geo- 
metrical distance between the geometrical points 2 and 1 
on a radial geodesic line in the frame of free falling clocks 
0- Nevertheless, even in this frame the absolute value 
of the variable p remains not fixed by the 3D distance 
measurements. 

In the static spherically symmetric case the choice 
of spherical coordinates and static metric dictates the 
form of three of the gauge fixing coefficients J2J): f t = 



o, r £ 



/Sit P 

are still not fixed 



cot#, = 0, but the form of the quantity 
and, equivalently, the function p(r) 

Here and further on, the prime de- 
notes differentiation with respect to the variable r. We 
refer to the freedom of choice of the function p(r) as " a 
rho-gauge freedom" , and to the choice of the p(r) function 
as " a rho-gauge fixing" in a broad sense. 

In the present article we will not use more general 
gauge transformations in a broad sense, than the rho- 
gauge ones. In our ID approach to the problem at hand 
all possible mathematical complications, due to the use 
of such wide class of transformations, can be easily con- 
trolled. 

The overall action for the aggregate of a point particle 
and its gravitational field is Atot = Aqr + Am - Ne- 
glecting the surface terms one can represent the Hilbert- 
Einstein action Agr — — 16 ^ G J d 4 Xy/\g\R and the me- 
chanical action Am = — -Wo J ds of the point source with 
a bare mass Mq as integrals with respect to the time 
t and the radial variable r of the following Lagrangian 
densities: 



£-gr. : 

£m 



f 



W (y/B)' '+{P') 2 V9Tt 



2Gn \ V~9r 

= -M ^gT t S(r). 



(4) 



Here Gm is the Newton gravitational constant, S(r) is 
the ID Dirac function (We are using units c=l.) 

As a result of the rho-gauge freedom the field vari- 
able \J—g rr is not a true dynamical variable but rather 



plays the role of a Lagrange multiplier, which is needed 
in a description of a constrained dynamics. This auxil- 
iary variable enters the Lagrangian Cgr in a nonlinear 
manner, in a contrast to the case of the standard La- 
grange multipliers. The corresponding Euler-Lagrange 
equations read: 
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9ttV-9rr = (5) 



where the symbol " = " denotes a weak equality in the 
sense of the theory of constrained dynamical systems. 

If one ignores the point source of the gravitational field, 
thus considering only the domain r > where S(r) = 0, 
one obtains the standard solution of this system Q : 

9tt(r) = 1- Po/p{r), g tt (r) g rr {r) = - (p(r)') , (6) 

where pa = 2GnM is the Schwarzschild radius, M is the 
gravitational (Keplerian) mass of the source, and p(r) is 
an arbitrary C 1 function. 



III. SOME EXAMPLES OF DIFFERENT 
RADIAL GAUGES 

In the literature one can find different choices of the 
function p(r) for the problem at hand: 

1. Schwarzschild gauge Q : p(r) = (r z + p'q) 1 ^ . It pro- 
duces f r = -2 (/o(r- 3 -^)-p G (r 3 /2-^))/(p 3 (/>-pG)). 

2. Hubert gauge [T(|: p(r) — r. It gives T r = 

~r i 1 ~ 2?) / i 1 ~ This sim P le choice of thc fac- 
tion p(r) is often related incorrectly with the original 

Schwarzschild article [5( (See [Tj])- In this case the co- 
ordinate r coincides with the luminosity distance p and 
the physical domain r € [0, oo) contains an event horizon 
at p H — pa- This unusual circumstance forces one to de- 
velop a nontrivial theory of black holes for Hilbert gauge 
- see for example 0, Il2l and the references therein. 

3. Droste gauge |13|: The function p(r) is given 
implicitly by the relations p/pc — cosh 2 -0 > 1 and 
t Ipg — ip + sinh^ cosh-0. The coefficient r. r = 

For this solution the vari- 



_± IP 1PG_ \ I /i _ £o 

a \ pa 4 I 1 \l 



4 p 



able r has a clear geometrical meaning: it measures the 
3D-radial distance to the center of the spherical symme- 
try. 

4 



Weyl gauge Q: p(r) = \ Ur] pc + Vpa/r) 
p G and f r = -f/(l- 4 



> 



__L/(^i_££y i n this 

PG 1 \PG r J 

gauge the 3D-space -M^ {~g mn (r)} becomes obviously 
conformally flat. The coordinate r is the radial variable 
in the corresponding Euclidean 3D space. 
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5. Einstein- Rosen gauge In the original article 
[H| the variable u 2 = p — pa has been used. To have a 
proper dimension we replace it by r — u 2 > 0. Hence, 
P(r) =r + PG >p G and f r = -f r -±^. 

6. Isotropic t-r gauge, defined according to the for- 
mula r = p+pa hi — lj , p > pa- Then only the com- 
bination (dt 2 — dr 2 ) appears in the 4D-interval ds 2 and 

7. One more rho-gauge in GR was introduced by Pu- 
gachev and Gun'ko @ and, independently, by Menzel 
[l6j. In comparison with the previous ones it is simple 
and more natural with respect to the quantity T r . We fix 
it using the condition f r — — - which is identical with the 
rho-gauge fixing for spherical coordinates in a flat space- 
time. Thus the curved space-time coordinates t, r, 6, (f> 
are fixed in a complete coherent way with the flat space- 
time spherical ones. 

Then, absorbing in the coordinates' units two inessen- 
tial constants, we obtain a novel form of the 4D interval 
for a point source: 

where N(r) = (2ip N )~ 1 (e 2 ^« - l); AT(r)~ l+0{^-) - for 
r^oo and N(r)~ - for In this specific gauge 

p(r) = r/N(r) and p(r) ~ r - for r — > oo, p(r) —> Pg ~ for 

The most remarkable property of this flat space-time- 
coherent gauge is the role of the exact classical New- 
ton gravitational potential ip N (r) = — ^ in the above 
exact GR solution. As usual, the component g u = 
e 2 ^»( r ) ~ 1 + 2tp N (r) + 0((p N (r)) - for r ^ oo. But 
now gtt has an essentially singular point at r = 0. 

The new form of the metric is asymptotically flat 
for r — > oo. For r — ► +0 it has a limit d0 2 + 
sin 2 8d(j) 2 . The geometry is regular in the whole space- 
time M ( ' 1, ' i \g tiu {x)) because the two nonzero scalars 

±R^x k R^ Xk = i(l-Soo) 4 - i(l-e 2 ^«) 4 and 

^R^X K R XK arR aT ^ = 4 (l-500) 6 = -| (l-e 2 ^«) 6 

are finite everywhere, including the center r = +0. At 
this center we have a zero 4D volume, because \f\g~\ = 
e 4v « ( r )jV(r) _8 r 2 sin 6*. At the same time the coordinates 
x m € (—00,00), to = 1,2,3; are globally defined in 
M {1 ' 3) (g^(x)). One can check that A g x m = 2 (f f x m . 
Hence, A g x m - O (^) - for r -> 00 and A g x m - 2x m ~ 
O (r) - for r — > +0, i.e., the coherent coordinates a;™ 1 are 
asymptotically harmonic in both limits, but not for finite 
values of r w pc- 

At this point one has to analyze two apparent facts: 
i) An event horizon p H exists in the physical domain 
only under Hilbert choice of the function p(r) = r, but 
not in the other gauges, discussed above. This demon- 
strates that the existence of black holes strongly depends 
on this choice of the rho-gauge in a broad sense. 



ii) The choice of the function p(r) can change drasti- 
cally the character of the singularity at the place of the 
point source of the metric field in GR. 

New interesting gauge conditions for the radial vari- 
able r and corresponding gau ge t ransformations were in- 
troduced and investigated in In the recent articles 
|18| it was shown that using a nonstandard p-gauge for 
the applications of GR to the stelar physics one can de- 
scribe extreme objects with arbitrary large mass, density 
and size. In particular, one is able to shift the value of 
the Openheimer-Volkoff limiting mass 0.7M Q of neutron 
stars to a new one: 3.8M©. This may shed a new light on 
these astrophysical problems and once more shows that 
the choice of rho-gauge (i.e. the choice of the coordinate 
r) can have real physical consequences. 

Thus we see that by analogy with the classical elec- 
trodynamics and non-abelian gauge theories of general 
type, in GR we must use indeed a more refined terminol- 
ogy and corresponding mathematical constructions. We 
already call the choice of the function p(r) "a rho-gauge 
fixing in a broad sense" . Now we see that the different 
functions p(r) may describe different physical solutions 
of Einstein equations Q with the same spherical sym- 
metry in presence of only one singular point at the center 
of symmetry. As we have seen, the mathematical prop- 
erties of the singular point may depend on the choice of 
the rho-gauge in a broad sense. Now the problem is to 
clarify the physical meaning of the singular points with 
different mathematical characteristics. 

Some of the static spherically symmetric solutions 
with a singularity at the center can be related by reg- 
ular gauge transformations. For example, choosing the 
Hilbert gauge we obtain a black hole solution of the vac- 
uum Einstein equations . Then we can perform a reg- 
ular gauge transformation to harmonic coordinates, de- 
fined by the relations T t — 0, r x i = 0, T x 2 = 0, r x 3 = 0, 
preserving the existence of the event horizon, i.e. staying 
in the same gauge sector. In the next sections we will 
find a geometrical criteria for answering the question: 
"When the different representations of static spherically 
symmetric solutions with point singularity describe dif- 
feomorphic space-times ?" 

The definition of the regular gauge transformations al- 
lows transformations that do not change the number and 
the character of singular points of the solution in the real 
domain and the behavior of the solution at the bound- 
ary, but can place these points and the boundary at new 
positions. 



IV. GENERAL FORM OF THE FIELD 
EQUATIONS IN HILBERT GAUGE AND SOME 
BASIC PROPERTIES OF THEIR SOLUTIONS 

According to Lichnerowicz 

IH E3, the physical 
space-times M.^ 1,3 '{g^, v {x)} of general type must be 
smooth manifolds at least of class C 2 and the metric 
coefficients g^ix) have to be at least of class C 3 , i.e. at 
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least three times continuously differentiable functions of 
the coordinates x. When one considers gtt{ r ), g rr (r) and 
p(r) not as distributions, but as usual functions of this 
class, one is allowed to use the rules of the analysis of 
functions of one variable r. Especially, one can multiply 
these functions and their derivatives of the correspond- 
ing order, one can raise them at various powers, define 
functions like log, exp and other mathematical functions 
of gttir), g r r(f), p(r) and coresponding derivatives. In 
general, these operations are forbidden for distributions 



In addition, making use of the standard rules for op- 
erating with Dirac 5-function Q and accepting (for sim- 
plicity) the following assumptions: 1) the function p(r) is 
a monotonic function of the real variable r; 2) p(0) = pa 
and the equation p(r) — po has only one real solution: 
r = 0; one obtains the following alternative form of the 
Eq.©: 



-In p 



1 



d 2 In g u , 1 fdlngtt 



(dlnp) 2 



dlnp 



d\np \ \-gpp 

2a G N M 5(p-p ), 
1 dlng tt \ d\ng pp 



2 d\np 
d\ng u 
dhip 



d\np 
9pp + 1 



= 0, 
= 0. (7) 



Note that the only remnant of the function p(r) in the 
system Eq.Q are the numbers po and oo = sign(p'(0)), 
which enter only the first equation, related to the source 
of gravity. Here gu and g pp — g rr / (p') 2 are considered as 
functions of the independent variable p G [po,oo). 

The form of the Eq. J7J shows why one is tempted to 
consider the Hilbert gauge as a preferable one: In this 
form the arbitrary function p(r) " disappears" . 

One usually ignores the general case of an arbitrary 
value po =/= accepting, without any justification, the 
value po = 0, which seems to be natural in Hilbert gauge. 
Indeed, if we consider the luminosity distance as a mea- 
sure of the real geometrical distance to the point source 
of gravity in the 3D space, we have to accept the value 
po = for the position of the point source. Otherwise 
the (5-function term in the Eq. (JJJ will describe a shell 
with radius po ^ 0, instead of a point source. 

Actually, according to our previous consideration the 
point source has to be described using the function 5(r). 
The physical source of gravity is placed at the point r = 
by definition. There is no reason to change this initial po- 
sition of the source, or the interpretation of the variables 
in the problem at hand. To what value of the luminosity 
distance po = p(0) corresponds the real position of the 
point source is not known a priori. This depends strongly 
on the choice of the rho-gauge. One can not exclude such 
a nonstandard behavior of the physically reasonable rho- 
gauge function p{r) which leads to some value po ^ 0. 
Physically this means that instead to infinity, the lumi- 
nosity of the point source will go to a finite value, when 



the distance to the source goes to zero. This very in- 
teresting new possibility appears in curved space-times 
due to their unusual geometrical properties. It may have 
a great impact for physics and deserves further careful 
investigation. 

In contrast, if one accepts the value po — 0, one has to 
note that the Hilbert-gauge singularity at p = will be 
space-like, not time-like. This is a quite unusual nonphys- 
ical property for a physical source of any static physical 
field. 

The choice of the values of the parameters po and o~o 
was discussed in |22j. There, a new argument in favor 
of the choice po — Pg and uq = 1 was given for differ- 
ent physical problems with spherical symmetry using an 
analogy with the Newton theory of gravity. In contrast, 
in the present article we shall analyze the physical mean- 
ing of the different values of the parameter po > 0. They 
turn out to be related with the gravitational mass defect. 

The solution of the subsystem formed by the last two 
equations of the system Q is given by the well known 
functions 



gtt(p) = l-pG/p, 9 P p(p)=-l/9tt{p)- 



(8) 



Note that in this subsystem one of the equations is a 
field equation, but the other one is a constraint. However, 
these functions do not solve the first of the Eq. (JJJ for any 
value of po, if Mq ^ 0. Indeed, for these functions the 
left hand side of the first field equation equals identically 
zero and does not have a 5-function-type of singularity, 
in contrast to the right hand side. Hence, the first field 
equation remains unsolved. Thus we see that: 

1) Outside the singular point, i.e. for p > po, the 
Birkhoff theorem is strictly valid and we have a stan- 
dard form of the solution, when expressed through the 
luminosity distance p. 

2) The assumption that <?«(?*) 5 grr{r) and p(r) are 
usual C 3 smooth functions, instead of distributions, 
yields a contradiction, if Mo ^ 0. This way one is not 
able to describe correctly the gravitational field of a mas- 
sive point source of gravity in GR. For this purpose the 
first derivative with respect to the variable r at least of 
one of the metric coefficients g^ v must have a finite jump, 
needed to reproduce the Dirac 5-function in the energy- 
momentum tensor of the massive point particle. 

3) The widespread form of the Schwarzschild's solution 
in Hilbert gauge JHJ does not describe a gravitational 
field of a massive point source and corresponds to the 
case Mq = 0. It does not belong to any physical gauge 
sector of gravitational field, created by a massive point 
source in GR. Obviously, this solution has a geometrical- 
topological nature and may be used in the attempts to 
reach pure geometrical description of "matter without 
matter" . In the standard approach to this solution no 
bare mass distribution was ever introduced. 

The above consideration confirms the conclusion, 
which was reached in |19| using isotropic gauge for 
Hilbert solution. In contrast to [lJj , in the next Sections 
we will show how one can include in GR neutral particles 
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with nonzero bare mass Mq. The corresponding new so- 
lutions are related to the Hilbert one via singular gauge 
transformations. 

V. NORMAL COORDINATES FOR 
GRAVITATIONAL FIELD OF A MASSIVE 
POINT PARTICLE IN GENERAL RELATIVITY 

An obstacle for the description of the gravitational field 
of a point source at the initial stage of development of 
GR was the absence of an adequate mathematical formal- 
ism. Even after the development of the correct theory of 
mathematical distributions Q there still exist an opin- 
ion that this theory is inapplicable to GR because of the 
nonlinear character of Einstein equations (|T|l . For exam- 
ple, the author of the article [21| emphasizes that "the 
Einstein equations, being non-linear, are defined essen- 
tially, only within framework of functions. The function- 
als, introduced in ... physics and mathematics (Dirac's 
5- function, "weak" solutions of partial differential equa- 
tions, distributions of Scwartz) are suitable only for linear 
problems, since their product is not, in general, defined." 
In the more recent article [2^ the authors have consid- 
ered singular lines and surfaces, using mathematical dis- 
tributions. They have stressed, that "there is apparently 
no viable treatment of point particles as concentrated 
sources in GR". See also Q and the references therein. 
Here we propose a novel approach to this problem, based 
on a specific choice of the field variables. 

Let us represent the metric ds 2 (QJ) of the problem at 
hand in a specific form: 

e 2vi dt 2 -e- 2vi+4v2 - 2 ^dr 2 -fe- 2vi+2v2 {d9 2 +sm 2 9d(/) 2 )(9) 

where ipi (r), <f2{f) and (p(r) are unknown functions of the 
variable r and p is a constant - the unit for luminosity 
distance p — p e~ vl+V2 . By ignoring the surface terms in 
the corresponding integrals, one obtains the gravitational 
and the mechanical actions in the form: 

A GR=2^J dt J &{f HM) 2 + (M) 2 )+^e 2v2 ), 
A Mo = - JdtjdrM e Vl S(r). (10) 

Thus we see that the field variables fi{r), fiir) and 
<p(r) play the role of a normal fields' coordinates in our 
problem. The field equations read: 

P 2 

A^ 2 (r) = i^2(r)-^)) (n) 
P 

where A r = e~^-^ (e^-^) is related to the radial part 
of the 3D-Laplacean, A = -l/ VPtflft (VWVl^dj) = 
-g» - f i d j )=A r + l/p 2 A e4l : A r = -g rr A r . 



The variation of the total action with respect to the 
auxiliary variable (p gives the constraint: 

(-(M) 2 + (M) 2 ) " H 0. (12) 

One can have some doubts about the correctness of the 
above derivation of the field equations l|ll|) , because here 
we use the Weyl's trick, applying the spherical symme- 
try directly to the action functional, not to the Einstein 
equations QJ. The correctness of the result of this pro- 
cedure is proved in the Appendix A. Therefore, if one 
prefers, one can consider the Lagrangian densities |@J, 
or the actions 1|10[) as an auxiliary tools for formulating 
of our one-dimensional problem, defined by the reduced 
spherically symmetric Einstein equations Q, as a vari- 
ational problem. The variational approach makes trans- 
parent the role, in the sense of the theory of constrained 
dynamical systems, of the various differential equations, 
which govern the problem. 

VI. REGULAR GAUGES AND GENERAL 
REGULAR SOLUTIONS OF THE PROBLEM 

The advantage of the above normal fields' coordinates 
is that when expressed through them the field equations 
(fTTf) are linear with respect to the derivatives of the un- 
known functions ^(r). This circumstance legitimates 
the correct application of the mathematical theory of dis- 
tributions and makes our normal coordinates privileged 
field variables. 

The choice of the function (p(r) fixes the rho-gauge in 
the normal coordinates. We have to choose this func- 
tion in a way that makes the first of the equations l|ll|) 
meaningful. Note that this non-homogeneous equation is 
quasi-linear and has a correct mathematical meaning if, 
and only if, the condition |y>i(0) — </?(0)| < 00 is satisfied. 

Let's consider once more the domain r > 0. In 
this domain the first of the equations (|11|) gives 
<pi{r) — C\ j e^^ r 'dr + Ci with arbitrary constants 
C± t 2' Suppose that the function (p(r) has an asymp- 
totics exp(— (p{r)) ~ kr n in the limit r — > +0 (with some 
arbitrary constants k and n). Then one easily obtains 
</?i(r) - (p{r) ~ Cikr n+1 /(n + 1) + ralnr + Infe + C 2 - if 
n ^ -1, and ipi(r) - (p(r) ~ (Cifc - 1) lnr + lnfc + C 2 - 
for n = — 1. Now we see that one can satisfy the condi- 
tion lim r ^o lvi( r ) — ¥?( r )l = constant < 00 for arbitrary 
values of the constants G\$ if, and only if n = 0. This 
means that we must have (p(r) ~ k — const ^ ±00 for 
r — > 0. We call such gauges regular gauges for the prob- 
lem at hand. Then <^i(0) = const ^ ±00. Obviously, the 
simplest choice of a regular gauge is (p(r) = 0. Further 
on we shall use this basic regular gauge. Other regu- 
lar gauges for the same gauge sector defer from it by a 
regular rho-gauge transformation which describes a dif- 
feomorphism of the fixed by the basic regular gauge man- 
ifold A4^{g mn (r)}. In terms of the metric components 
the basic regular gauge condition reads p 4 gtt + p 4 grr = 
and gives r r = 0. 
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Under this gauge the field equations (|llfl acquire the 
simplest quasi-linear form: 

G^Mo^ i(0) L e ^(r)_ (13) 

P P 

The constraint (|f 2() acquires the simple form: 



(M) 2 + (pV 2 ) 2 -e 2 ^=0. 



(14) 



As can be easily seen, the basic regular gauge (p(r) = 
has the unique property to split the system of field equa- 
tions (|ll|l and the constraint l|12() into three independent 
relations. 

The new field equations ((13(1 have a general solution 



ip 2 (r) = -\n 



1 . / ^ — roo — r 
sinn V2e 



(15) 



The first expression in Ea. (|15l) represents a distribution 
tfi(r). In it O(r) is the Heaviside step function. Here we 
use the additional assumption 8(0) := 1. It gives a spe- 
cific regularization of the products, degrees and functions 
of the distribution 0(r) and makes them definite. For ex- 
ample: (6(r)) 2 = 6(r), (6(r)) 3 = 6(r), . . . , Q(r)S(r) = 
5(r), _/(re(r))=/(r)e(r)-/(0)(6(r)-e(0))- for any 
function /(r) with a convergent Taylor series expansion 
around the point r = 0, and so on. This is the only sim- 
ple regularization of distributions we need in the present 
article. 

The second expression <p2(r) in Ea. (|15|) is a usual 
function of the variable r. The symbol is used 
as an abbreviation for the constant expression = 
sign (</? 2 (0)) p sinh (^/Si^e" 1 ^ )) /V^2- 

The constants 

ei = -^- 2 ^iM 2 +^^(o)e^(°)(e(r)-e(o)), 



are the values of the corresponding first integrals l|l(j|) of 
the differential equations (|13fl for a given solution (|15l) . 

Then for the regular solutions (|15f) the condition 114|) 
reads: 



(e(r)-e(o)) = o. (i7) 



An unexpected property of this relation is that it can- 
not be satisfied for any value of the variable r £ (— oo, oo), 
because £1^2 are constants. The constraint (|17|) can be 
satisfied either on the interval r £ [0,oo), or on the in- 
terval r £ (—oo,0). If, from physical reasons we chose it 
to be valid at only one point r* £ [0,oo), this relation 
will be satisfied on the whole interval r £ [0, oo) and this 
interval will be the physically admissible real domain of 
the radial variable. Thus one can see that our approach 



gives a unique possibility to derive the admissible real 
domain of the variable r from the dynamical constraint 
(fpf)l . i.e., this dynamical constraint yields a geometrical 
constraint on the values of the radial variable. 

As a result, in the physical domain the values of the 
first integrals ((16(1 are related by the standard equation 



£tot — £i + £2 



0. 



(18) 



which reflects the fact that our variation problem is in- 
variant under re-parametrization of the independent vari- 
able r. At the end, as a direct consequence of the rela- 
tion (|18|) one obtains the inequality £2 — —£% > 0, be- 
cause in the real physical domain r £ [0,oo) we have 



^p 2 (Pi(r) = const < 0. 



For the function p r eg,o( r ) > 0, which corresponds to 
the basic regular gauge, one easily obtains 

-1 



Preg\e =0 {r) = PG 1 ~ exp 



(^)) 



(19) 



Now one has to impose several additional conditions 
on the solutions 1(15(1 : 

i) The requirement to have an asymptotically flat 
space-time. The limit r ~ » r^ corresponds to the 
limit p — > 00. For solutions l|15J) we have the property 
9rr{f 00) / p' {foa) = 1- The only nontrivial asymptotic 
condition is gtt{r oo) = 1. It gives tp' 1 {0)r oo -\-tpi{0) = 0. 

ii) The requirement to have the correct Keplerian mass 
M, as seen by a distant observer. Excluding the variable 
r > from g tt (r) = e 2 ^'') and p(r) = pe-^M+^M for 
solutions HI 5(1 one obtains gu = 1 + co " st , where const = 
2sign{p) sign(<p' 2 (0)) p 2 <^i(0) = -2G N M. 

iii) The consistency of the previous conditions with 
the relation g u g T r + p' 2 = gives signfroo) = sign(p) = 
sign(<p[(0)) = sign{cp' 2 (G}) = 1. 

iv) The most suitable choice of the unit p is p = 
GjyM = pg/2- As a result all initial constants become 
functions of the two parameters in the problem - 

and M: ^(0) = p 2 (0) = -m(sinh(^ 



(16) ^(0) 



n ^ a (0) 



■ coth 



G N M' r 2 \"J — G N M 

v) The gravitational defect of the mass of a point par- 
ticle. 

Representing the bare mechanical mass Mq of the 
point source in the form M — J^°° M S(r)dr = 

4tt J r °° y / —g rr (r)p 2 (r)iJ,(r)dr, one obtains for the 
mass distribution of the point particle the expression 

p(r) = M 5(r)l (47r A /- ffrr (r) p 2 (r)) = M S g {r), where 

S g (r) := 5(r)/ (^in \J — g rr (r) p 2 (r)j is the ID invariant 

Dirac delta function. The Keplerian gravitational mass 
M can be calculated using the Tolman formula yj: 



M = 47r/ pf(r)/?(r)ii(r)dr = M o ^g tt (0). (20) 
Jo 



Here we use the relation p 1 = y/—g~tt 9rr- As a result we 
reach the relations: gu(0) = e 



= p 2 Vl (Q) 



I ' x : ) 



( 2 G T nm) 
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< 1 and r x = G N M]n{^f) > 0. (Note that 

due to our convention 9(0) := 1 the component gtt(r) 
is a continuous function in the interval r G [0, oo) and 
5tt(0) = 5tt(+0) is a well defined quantity.) 

The ratio g = = \fg~ttW) G [0, 1] describes the 
gravitational mass defect of the point particle as a second 
physical parameter in the problem. The Keplerian mass 
M and the ratio g define completely the solutions (|15|l . 

For the initial constants of the problem one obtains: 



ifi(0) = lng, tp 2 (Q) = - In 
1 



i-e 2 



G N M 



2g ' 
1 1 + g 2 



G N M l-g 2 ' 



(21) 



Thus we arrive at the following form of the solutions 
G3: 



r Q(r) 



Pa (r) = - In -(]/<? 



G N M 

,-r/G N M _ np r/G N M 



g e 



f )) ( 22 ) 



and the rho-gauge fixing function 



Preg\e= (r) = Pg (1 - Q exp ( 4r/p G 



(23) 



An unexpected feature of this £wo parametric variety 
of solutions for the gravitational field of a point particle is 
that each solution must be considered only in the domain 
r G [0, GnMIii(1/ g)], if we wish to have a monotonic 
increase of the luminosity distance in the interval [po, oo). 

It is easy to check that away from the source, i.e., for 
r > 0, these solutions coincide with the solution JHJ in 
the Hilbert gauge. Hence, outside the source the solu- 
tions l|22l) acquire the well known standard form, when 
represented using the variable p. This means that the 
solutions l|22|) strictly respect a generalized Birkhoff the- 
orem. Its proper generalization requires only a justifi- 
cation of the physical domain of the variable p. It is 
remarkable that for the solutions (|22|l the minimal value 
of the luminosity distance is 



po = 2G N M/(1 - g 2 ) > PG . 



(24) 



This changes the Gauss theorem and leads to important 
physical consequences. One of them is that one must ap- 
ply the Birkhoff theorem only in the interval p G [po, oo). 



VII. REGULAR GAUGE MAPPING OF THE 
INTERVAL r G [0, r<x>] ONTO THE WHOLE 
INTERVAL r G [0, oo] 

As we have stressed in the previous section, the so- 
lutions 11' I'D for the gravitational field of a point par- 
ticle must be considered only in the physical domain 



r G [0, loo]. It does not seem to be very convenient to 
work with such unusual radial variable r. One can eas- 
ily overcome this problem using the regular rho-gauge 
transformation 



r/f 
<•/? + 1 



(25) 



with an arbitrary scale f of the new radial variable r 
(Note that in the present article we are using the same 
notation r for different radial variables.) This linear frac- 
tional diffeomorfism does not change the number and the 
character of the singular points of the solutions in the 
whole compactified complex plane Cr of the variable r. 
The transformation l|25[) simply places the point r = 
at the infinity r — oo, at the same time preserving the 
initial place of the origin r = 0. Now the new variable 
r varies in the standard interval r G [0, oo), the regular 
solutions 11' I'D acquire the final form 

Mr) = HVq) (i - ^TT e ( r / f )) ' 

Mr) = - In Q ((l/g) 7 ^ - g^hr^ , 

<p(r) = 2 ln(r/r + 1) + Xn^/r^). (26) 
The final form of the rho-gauge fixing function reads: 



Preg(r) = /0(J ( 1 - g^Tl 



(27) 



The last expression shows that the mathematically ad- 
missible interval of the values of the ratio g is the open 
interval (0,1). This is so, because for g = and for 
g = 1 we would have impermissible trivial gauge-fixing 
functions p reg (r) = l and p reg (r) = 0, respectively. 

Now we are ready to describe the singular character 
of the coordinate transition from the Hilbert form of 
Schwarzschild solution © to the regular one (|26fl . To 
simplify notations, let us introduce dimcnsionlcss vari- 
ables z = p/p g and C — r/f. Then the Eq. I|27|l shows 
that the essential part of the change of the coordinates 
is described, in both directions, by the functions: 



C(z) = 



40 = - 

ln(l/g 2 ) 
In z — ln(z — 1 



0< 



and 



-1, 6 6(0,1). 



(28) 



Obviously, the function z(() is regular at the place of 
the point source £ = 0; it has a simple pole at £ = oo and 
an essential singular point at £ = — 1. At the same time 
the inverse function £(z) has a logarithmic branch points 
at the Hilbert-gauge center of symmetry z = and at the 
"event horizon" 2 = 1. Thus we see how one produces 
the Hilbert-gauge singularities at p = and at p = p H , 
starting from a regular solution. The derivative 



dC/dz 



ln(l/g 2 ) 



z[z- l)(lnz-ln(z- 1))" 
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aproaches infinity at these two points, hence the singular 
character of the change in the whole complex domain 
of the variables. Thus we reach a complete description 
of the change of coordinates and its singularities in the 
complex domain of the radial variable. 

The restriction of the change of the radial variables on 
the corresponding physical interval outside the source: 
C G (0, oo) hz£ (1/(1 — g 2 ), oo), is a regular one. 

The expressions iffifl) and lj27|) still depend on the 
choice of the units for the new variable r. We have 
to fix the arbitrary scale of this variable in the form 
f = p G / ln(l/g 2 ) = G N M/ In (^) to ensure the validity 
of the standard asymptotic expansion: g^ ~ 1 — pq fr + 
O ((pg/ 1 ") 2 ) when r — > oo. 

Then the final form of the 4D interval, defined by the 
new regular solutions outside the source is: 



ds 2 = e 2 *a dt 2 



dr 2 



^-^J-Krf (^ 2 + sinW). (29) 



Here we are using a modified (Newton-like) gravitational 
potential: 



<p G (r;M,M ) 



a coefficient N G (r) — (2<^ G )~ 
luminosity distance 



G N M 



G N M/\n(Mz 



(30) 



l) , and an optical 



, s r+G N M/\n(^) 
p[r) =2G N M/(l- e 2 *o ) = ' 



(31) 



These basic formulas describe in a more usual way our 
regular solutions of Einstein equations for r S (0, oo). 
They show immediately that in the limit g — > our so- 
lutions tend to the Pugachev-Gun'ko-Menzel one, and 
ip G (0;M, Mo) = lag — ► — oo. In the limit: g — > 1 we 
obtain for any value of the ratio r/p g : gtt(r j pa , g) — * 1, 
g rr (r/pG, g) — * — 1 5 and p(r/pc, g) — * oo. Because of the 
last result the 4D geometry does not have a meaningful 
limit when g — > 1. In this case <p G (r;M,Mo) — » at 
all 3D space points. Physically this means that solutions 
without mass defect are not admissible in GR. 



VIII. TOTAL ENERGY OF A POINT SOURCE 
AND ITS GRAVITATIONAL FIELD 

In the problem at hand we have an extreme example of 
an "island universe". In it a privileged reference system 
and a well defined global time exist. It is well known 
that under these conditions the energy of the gravita- 
tional field can be defined unambiguously . Moreover, 
we can calculate the total energy of the aggregate of a 
mechanical particle and its gravitational field. 

Indeed, the canonical procedure produces a total 



Hamilton density Tttot 

1 ! -2/2,-2/2 



M e Vl S(r). Using the 



constraint (|18fl and the first of the relations (|21(l , one im- 
mediately obtains for the total energy of the GR universe 
with one point particle in it: 



Htotdr — M = qM < M . 



(32) 



This result completely agrees with the strong equiv- 
alence principle of GR. The energy of the gravitational 
field, created by a point particle is the negative quantity: 
E GR = Ekat -Eo = M - M = -Af (l -g)<0. 

The above consideration gives a clear physical explana- 
tion of the gravitational mass defect for a point particle. 



IX. INVARIANT CHARACTERISTICS OF THE 
SOLUTIONS WITH A POINT SOURCE 

A. Local Singularities of Point Sources 

Using the invariant ID Dirac function one can write 
down the first of the Eq. in a form of an exact 

relativistic Poisson equation: 



A r <pi(r) =4nG N M Q Sg{r). 



(33) 



This equation is a specific realization of Fock's idea (see 
Fock in U) using our normal fields coordinates (Sec. V). 
It can be re-written, too, in a transparent 3D form: 



Api(r) =47rGjvMo^ 3) (r 



(34) 



The use of the invariant Dirac 5 ff -function has the 
advantage that under diffeomorphisms of the space 
5mn(r)} the singularities of the right hand side of 
the relativistic Poisson equation l|33ll remain unchanged. 
Hence, we have the same singularity at the place of the 
source for the whole class of physically equivalent rho- 
gauges. Then one can distinguish the physically different 
solutions of Einstein equations with spherical symme- 
try by investigating the asymptotics in the limit r —>■ +0 

of the coefficient j(r) — 1/ ^47rp(r) 2 yj — g rr {r)j in front 

of the usual ID Dirac <5(r)-function in the representation 
5 g (r) = j(r)8(r) of the invariant one. 

For regular solutions (J2SJ| the limit r — > +0 of 

this coefficient is a constant 



7(0) 



1 



Airp z G \g 



- - <J 



(35) 



which describes the intensity of the invariant (5-function 
and leads to the formula 8 g (r) = y(0)6(r). We see that: 

1) The condition g £ (0, 1) ensures the correct sign of 
the intensity 7(0), i.e., the property 7(0) € (0,oo), and 
thus a negativity of the total energy Eqr of the gravi- 
tational field of a point particle, according the previous 
Section III. 

2) For g — we have 7(0) = oo, and for g=l: 7(0) = 0. 
Hence, we have non-physical values of the intensity 7(0) 
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in these two cases. The conclusion is that the physical 
interval of values of the ratio g is the open interval (0, 1). 
This is consistent with the mathematical analysis of this 
problem, given in the previous Section. 

It is easy to obtain the asymptotics of the coefficient 
7(r) in the limit r — > for other solutions, considered in 
the Introduction: 



Schwarzschild solution: 

Hilbert solution: 
Droste solution: 

Weyl solution: 
Einstein-Rosen solution: 

Isotropic (t-r) solution: 
Pugachev-Gun'ko-Menzel 
solution: 



7 S M ~ 

Id ( r ) ~ 
Iw 0) r 

7 M (r) 

1er{t) 

IpGM V) ' 



1 (PG^ 1 ! 2 . 
4np' G V r ) ' 

i (pg_\ 5 / 2 . 
4irp G V r J ' 
1 . 



16 ( r \ 

1 (PG_ 

7rp G \ r 



1 (PG_\ 1 / 2 . 

4ttp g V r J ' 
1/2 



p G /2r 



As we can see: 

1) Most of the listed solutions are physically differ- 
ent. Only two of them: Schwarzschild and Einstein- 
Rosen ones, have the same singularity at the place of the 
point source and, as a result, have diffeomorphic spaces 
A^ 3 Mff«m( r )}> which can be related by a regular rho- 
gauge transformation. 

2) As a result of the alteration of the physical meaning 
of the variable p = r inside the sphere of radius p H , 
in Hilbert gauge the coefficient 7(7") tends to imaginary 
infinity for r — ► 0. This is in a sharp contrast to the real 
asymptotic of all other solutions in the limit r — > 0. 

3) The Droste solution is a regular one and cor- 
responds to a gravitational mass defect with g D — 
(v%-l)/2 fa 0.61803. The golden-ratio-conjugate 

number (\/5 — l) /2, called sometimes also "a silver ra- 
tio" , appears in our problem as a root of the equation 
1/ q—q = 1, which belongs to the interval g £ (0, 1). (The 
other root of the same equation is — (a/5 + l) /2 < and 
does not correspond to a physical solution.) Hence, the 
Droste solution can be transformed to the form l|22|l . or 
(|26[) by a proper regular rho-gauge transformation. 

4) The Weyl solution resembles a regular solution with 
g w = 1 and 7 W (0) = 0. Actually the exact invariant 
(5-function for this solution may be written in the form 
(see the Appendix B): 



J g,w 



(r) 
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This distribution acts as a zero-functional on the stan- 
dard class of smooth test functions with compact sup- 
port, which are finite at the points r = 0, 00. 
The formula (|36|) shows that: 

a) The Weyl solution describes a problem with a spher- 
ical symmetry in the presence of two point sources: one 
with grL = 1 at the point r = and another one with 



g^, = 1 at the point r = 00. Hence, it turns out that 
this solution is the first exact analytical two-particle-likc 
solution of Einstein equations JQ. Unfortunately, as a 
two-point solution, the Weyl solution is physically trivial: 
it does not describe the dynamics of two point particles 
at a finite distance, but rather gives only a static state of 
two particles at an infinite distance. 

b) To make these statements correct, one has to com- 
pactify the comformally flat 3D space A4^{— g mn (r)} 
joining it to the infinite point r = 00. Then both the 
Weyl solution and its source will be invariant with re- 
spect to the inversion r — ► p%jr. 

5) The isotropic (t-r) solution resembles a regular one, 
but in it x ps 1.27846 is the only real root of the equation 
x + hi(x — 1) = and gives a non-physical value of the 
parameter g T pa 2.14269, which does not belong to the 
physical interval (0, 1). 

6) The singularity of the coefficient 7(1") in front of 
the usual ID Dirac (5-function in the Pugachev-Gun'ko- 
Menzel solution is stronger than any polynomial ones in 
the other listed solutions. 



B. Local Geometrical Singularities of the Regular 
Static Spherically Symmetric Space-Times 

It is well known that the scalar invariants of the Rie- 
mann tensor allow a manifestly coordinate independent 
description of the geometry of space-time manifold. The 
problem of a single point particle can be considered as 
an extreme case of a perfect fluid, which consists of only 
one particle. According to the article j2^|, the maximal 
number of independent real invariants for such a fluid is 
9. These are the scalar curvature R and the standard 
invariants r\, u>i, W2, ms, 777,5. (See [24[ for their defini- 
tions.) The invariants i?, n and 7713 are real numbers 
and the invariants Wi,W2 and 777.5 are, in general, com- 
plex numbers. The form of these invariants for spheri- 
cally symmetric space-times in normal field variables is 
presented in Appendix C. Using these formulas we ob- 
tain for the regular solutions l|22l) in the basic gauge the 
following simple invariants: 



h = 



8P% Q 



-S{ — 

KPG 



1 (l-g 2 )\(r\ 7r(l-g 2 ) 3 



2 "V^; 2 g " 9 \p G J 2' 



1 



h = o, 

_ e(r/p G )-l f 2 r 4r/ p G \ 4 _ Pg e(r/ PG )-i 
3 8p 2 G g 2 \ J 8g 2 p(r) 4 ' 

©O/Pg) A g 2 c 4r/p G y _PG ®(r/p G ) 



4p 



p(rf 



(37) 



As can be seen easily: 

1) The invariants I\,...,a of the Riemann tensor are well 
defined distributions. This is in accordance with the gen- 
eral expectations, described in the articles [2!|, where 
one can find a correct mathematical treatment of distri- 
butional valued curvature tensors in GR. 
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As we saw, the manifold M^ 1 ^ {g^x)} for our regu- 
lar solutions has a geometrical singularity at the point, 
where the physical massive point particle is placed. 

Note that the metric for the regular solutions is glob- 
ally continuous, but its first derivative with respect to 
the radial variable r has a finite jump at the point source. 
This jump is needed for a correct mathematical descrip- 
tion of the delta- function distribution of matter of the 
massive point source in the right hand side of Einstein 
equations i ffl l. 

During the calculation of the expressions ^3,4 we have 
used once more our assumption 9(0) = 1. 

Three of the invariants Ij37|> are independent on the 
real axes r G (—00, 00) and this is the true number of the 
independent invariants in the problem of a single point 
source of gravity. On the real physical interval r £ [0, 00) 
one has I3 — and we remain with only two independent 
invariants. For r € (0, 00) the only independent invariant 
is I4, as is well known from the case of Hilbert solution. 

2) The scalar invariants ^i 4 have the same form as 
in Eq. (|37[1 for the regular solutions in the representation 
(|26|l . obtained via the diffeomorphic mapping i|25|) . 

3) All other geometrical invariants r 1 (r), Wi^(r), 
m 3 5 (r) include degrees of Dirac 5- function and are not 
well defined distributions. Therefore the choice of the 
simple invariants is essential and allows us to use 
solutions of Einstein equations, which are distributions. 

To the best of authors knowledge, the requirement to 
make possible the correct use of the theory of distribu- 
tions is a novel criteria for the choice of invariants of the 
curvature tensor and seems to not have been used until 
now. It is curious to know if it is possible to apply this 
new criteria in more general cases as opposed to the case 
of static spherically symmetric space-times. 

4) The geometry of the space-time depends essentially 
on both parameters M and g, which define the regular 
solutions (|26[l of Einstein equations Hence, for differ- 
ent values of the two parameters these solutions describe 
non-diffcomorphic space-times with different geometry. 



C. Event Horizon for Static Spherically Symmetric 
Space-Times with Point Singularity 

According to the well known theorems by Hawking, 
Penrose, Israel and many other investigators, the only 
solution with a regular event horizon not only in GR, 
but in the theories with scalar field(s) and in more gen- 
eral theories of gravity, is the Hilbert one 01 ■ As we 
have seen, this solution has a pure geometrical nature 
and does not describe a gravitational field of a point par- 
ticle with bare mechanical mass Mo 7^ 0. The strong 
mathematical results like the well known no hear theo- 
rems, etc, for the case of metrics with an event horizon 
can be shown to be based on to the assumption that such 
horizon is indeed present in the solution. These mathe- 
matical results do not contradict to the ones, obtained 
in the present article, because the other solutions, which 



we have considered together with the Hilbert one, do not 
have an event horizon at all. 

Indeed, for the point p H at which gtt(p H ) = one 
obtains p(p H ) = pa- The last equations do not have 
any solution p H 6 for the regular solutions i|22[l . 
(|26[) . The absence of a horizon in the physical domain 
r G [0,oo) is obvious in the representation (|29|) of the 
regular solutions. 

For the other classical solutions one obtains in the 
limit r — * p H as follows: 



Schwarzschild solution: 
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Einstein-Rosen solution: 
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As one can see, only for Hilbert, Schwarzschild and 
Einstein- Rosen solutions the metric component g r r{r) 
has a simple pole at the point p H . In the last two 
cases this is a point in the corresponding manifold 
Ai^ {(?mn(r)} , where the center of symmetry is placed, 
not a horizon. 



X. ON THE GEOMETRY OF MASSIVE POINT 
SOURCE IN GR 

It is curious that for a metric, given by a regular solu- 
tion, the 3D-volume of a ball with a small radius n, < pa, 
centered at the source, is 



Vol{r b ) 



3 (1 - e z Y p G p G 



It goes to zero linearly with respect to — > 0, in con- 
trast to the Euclidean case, where VolE(rb) ~ rf. This 
happens, because for the regular solutions l|22|) . (|26|l we 
; sin 9 — > po sin 9 7^ in the limit 



obtain -\/| 3 g| = p(r) 
r — > 0. Nevertheless, lim rb ->o Vol(rt) — and this legit- 
imates the use of the term " a point source of gravity" in 
the problem at hand: the source can be surrounded by a 
sphere with an arbitrary small volume Vol(rt) in it and 
with an arbitrary small radius r^. 

In contrast, when rb — > the area of the ball's surface 
has a finite limit: nz^p > 47rpg, and the radii of the 

big circles on this surface tend to a finite number ^% > 
2irp G . 

Such unusual geometry, created by the massive point 
sources in GR, may have an interesting physical con- 
sequences. For example, space-times, defined by the 
regular solutions lj2*2l . pp]l. have an unique property: 
When one approaches the point source, its luminosity 
remains finite. This leads to a very important modifica- 
tion of Gauss theorem in the corresponding 3D spaces. 
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After all, this modification may solve the well known 
problem of the classical divergences in field theory, be- 
cause, as we see, GR offers a natural cut-off parameter 
r = pa / ln(l/p 2 ) for fields, created by massive point par- 
ticles. If this will be confirmed by more detailed cal- 
culations, the price, one must pay for the possibility to 
overcome the old classical divergences problem will be to 
accept the idea, that the point objects can have a nonzero 
surface, having at the same time a zero volume and a zero 
size r. 

On the other hand, as we have seen, the inclusion of 
point particles in GR is impossible without such unusual 
geometry. The Hilbert solution does not offer any more 
attractive alternative for description of a source of grav- 
ity as a physical massive point, placed at the coordinate 
point p = 0, and with usual properties of a geometrical 
point. Hence, in GR we have no possibility to introduce 
a notion of a physical point particle with nonzero bare 
mass without a proper change of the standard definition 
for point particle as a mathematical object, which has 
simultaneously zero size, zero surface and zero volume. 
We are forced to remove at least one of these three usual 
requirements as a consequence of the concentration of the 
bare mass Mq at only one 3D space point. 

From physical point of view the new definition for 
physical point particle, at which we arrived in this ar- 
ticle, is obviously preferable. 

Moreover, it seems natural to have an essential differ- 
ence between the geometry of the "empty" space-time 
points and that of the matter points in GR. The same 
geometry for such essentially different objects is possi- 
ble only in the classical physics, where the space-time 
geometry does not depend on the matter. Only under 
the last assumption there is no geometrical difference be- 
tween pure mathematical space points and matter points. 
The intriguing new situation, described in this article, 
deserves further careful analysis. 

The appearance of the above nonstandard geometry in 
the point mass problem in GR was pointed at first by 
Marcel Brillouin |27| . Our consideration of the regular 
solutions (|22[l. (|26[) confirms his point of view on the 
character of the singularity of the gravitational field of 
massive point particles in GR. 



XI. CONCLUDING REMARKS 

The most important result of the present article is the 
explicit indication of the fact that there exist infinitely 
many different static solutions of Einstein equations 
with spherical symmetry, a point singularity, placed at 
the center of symmetry, and vacuum outside this sin- 
gularity, and with the same Keplerian mass M . These 
solutions fall into different gauge classes, which describe 
physically and geometrically different space-times. Some 
of them were discovered at the early stage of develop- 
ment of GR, but up to now they are often considered 
as equivalent representations of some "unique" solution 



which depends on only one parameter - the Keplerian 
mass M . As shown in Section IX A, this is not the case. 

As a consequence of Birkhoff theorem, when expressed 
in the same variables, the solutions with the same mass 
M indeed coincide in their common regular domain of 
coexistence - outside the singularities. In contrast, they 
may have a different behavior at the corresponding sin- 
gular points. 

A correct description of a massive point source of grav- 
ity is impossible, making use of most of these classical 
solutions. 

Using novel normal coordinates for the gravitational 
field of a single point particle with bare mechanical mass 
Mq we are able to describe correctly the massive point 
source of gravity in GR. The singular gauge transfor- 
mations yield the possibility to overcome the restriction 
to have a zero bare mass Mq for neutral point particles 
in GR ^ij. It turns out that this problem has a two- 
parametric family of regular solutions. 

One of the parameters - the bare mass Mq of the point 
source, can be obtained in a form of surface integral, 
integrating both sides of Eq. ffi^ on the whole 3D space 
M^{-g mn (r)}: 




The second parameter - the Keplerian mass M, is de- 
scribed, as shown in Section VIII, by the total energy: 

M= f d 3 rvfg~\Htot = [ d 3 r^\g~\ Mq 6® ( r ). 

Jm 3 Jm 3 

These two parameters define the gauge class of the 
given regular solution. Obviously the parameters M 
and Mq are invariant under regular static gauge trans- 
formations, i.e., under diffeomorphisms of the 3D space 
M^{— jmn(r)}. It is an analytical manifold with a 
strong singularity at the place of the massive point parti- 
cle. For every regular solution both parameters are finite 
and positive and satisfy the additional physical require- 
ment < M < Mq. 

It is convenient to use a more physical set of continuous 
parameters for fixing the regular solutions, namely: the 
Keplerian mass M £ (0,oo) and the gravitational mass 
defect ratio g = jf e (0, 1). 

The only classical solution, which is regular, is the 
Droste one. For this solution the gravitational mass de- 
fect ratio is g D — (\/5 — l) /2. 

For the regular solutions the physical values of the op- 
tical luminosity distance p are in the semi-constraint in- 
terval p G [j^,oo). 

Outside the source, i.e., for p > , the Birkhoff 
theorem is strictly respected for all regular solutions. 

The metric, defined by a regular solution, is globally 
continuous, but its first derivatives have a jump at the 
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point source. Our explicit results justify the Raju's con- 
clusion [2^] that GR will be consistent with the existence 
of point particles if one assumes the metric to be at most 
of class C° and show that the Hawking-Penrose singular- 
ity theory, based on the assumption that the components 
of g^v are at least C 1 functions, must be reconsidered. 

For the class of regular solutions, written in the form 
l)22[l . 1(26(1 . or l|29|l . the non-physical interval of the op- 
tical luminosity distance p € [0, j^s ], which includes 
the luminosity radius p H = pa, can be considered as an 
"optical illusion". All pure mathematical objects, which 
belong to this interval, are in the imaginary domain "be- 
hind" the real physical source of the gravitational field 
and have to be considered as a specific kind of optical 
"mirage". This is in agreement with Dirac's conclusion 
about the non-physical character of the inner domain 
P < Pg f° r the Hilbert solution [2(j and extends this 
conclusion on the whole non-physical interval of the op- 
tical luminosity distance p £ [0, jz^s] > P H f° r a given 
regular solution. 

We are forced to cut the Hilbert form of the 
Schwarzschild solution at the value po = jz^2 because 
of the presence of the matter point mass. Its presence ul- 
timately requires a definite jump of the first derivatives 
of the metric. The jump is needed to make the Einstein 
tensor coherent with the energy-momentum stress tensor 
of the point particle. This cutting can be considered as 
a further development of Dirac's idea pr| : 

"Each particle must have a finite size no smaller than 
the Schwarzschild radius. I tried for some time to work 
with a particle with radius equal to the Schwarzschild ra- 
dius, but I found great difficulties, because the field at 
the Schwarzschild radius is so strongly singular, and it 
seems that a more profitable line of investigation is to 
take a particle bigger than the Schwarzschild radius and 
to try to construct a theory for such particle interacting 
with gravitational field. " 

For a precise understanding of these statements one 
has to take into account that the above Dirac's idea is 
expressed in terms of the luminosity distance p. 

An alternative development, which is much closer to 
Dirac approach to the above problem, can be found in 
the recent articles pgj . 

The geometry around physical massive point parti- 
cles is essentially different from the geometry around the 
"empty" geometrical points. This unusual geometry, de- 
scribed at first by Brillouin, may offer a new way for 
overcoming of the classical fields divergence problems. 
The new possibility differs from the one, used in for 
charged point particles. It is based on the existence of a 
natural cut-off parameter f = G N M (c 2 In (^y-)) and 
on a new interpretation of the relation between the singu- 
lar mass distribution of a point particle and the geometry 
of the space-time around this particle. 

According to the remarkable comment by Poincare j2£j 
real problems can never be considered as solved or un- 
solved, but rather they are always more or less solved. 

The strong physical singularity at the "event horizon" 



of the Hilbert solution, stressed by Dirac and further 
physical consequences, which one can derive from the 
new regular solutions 1(26(1 of the old problem, consid- 
ered here will be discussed in a separate article. In this 
one we will add some more remarks. 

Our consideration shows that the observed by the as- 
tronomers compact dark objects (CDO), called by them 
black holes without any direct evidences for the existence 
of real event horizons [3J| ; can not be described theoreti- 
cally by solutions JHJ in Hilbert gauge, if we assume that 
these objects are made from matter of some nonzero bare 
mechanical mass Mo 7^ 0. At present the only real fact 
is the existence of a massive invisible CDO with Kep- 
lerian mass M, which is too large with respect to the 
conventional understanding of the stars physics. 

Concerning these unusual CDO, most probably one ac- 
tually has to solve a much more general problem. Since a 
significant amount of invisible dark matter, which mani- 
fests itself only due to its gravitational field, is observed 
in the Nature at very different scales: in the clusters of 
galaxies, in the halos of the galaxies, at the center of our 
galaxy, and as a compact dark components in some bi- 
nary star systems, one is tempted to look for some univer- 
sal explanation of all these phenomena. Obviously such 
universal explanation can not be based on Schwarzschild 
solution in Hilbert gauge and one has to look for some 
other theoretical approach. A similar idea was pointed 
out independently in the recent articles [3l| . 

At the end we wish to mention some additional open 
problems, both mathematical and physical ones, con- 
nected with the gravitational field of point particles. 

Our consideration was essentially restricted to the real 
domain of variables. The only exception was the de- 
scription of the singular change 1(28(1 of radial variable. 
It is well known that the natural domain for study 
of the solutions of holomorphic differential equations is 
the complex one. A complete knowledge of the solu- 
tions of such equations is impossible without descrip- 
tion of all singularities of the solutions in the complex 
domain. Therefore, looking for a complete analysis of 
the solutions of some differential equation of n-th order 
f (w^(z),w^-^{z),...,w^(z),w(z),z) = for a func- 
tion w(z),z S C^-\ where /(...) is a holomorphic func- 
tion of the corresponding complex variables, one has to 
consider these solutions as a holomorphic functions of 
complex variable z. Now the most important issue be- 
comes the study of the singular points of the solutions in 
the whole complex domain. 

For the Einstein's equations (JIJ such a four- 
dimensional complex analysis is impossible at present, 
because the corresponding mathematical methods are 
not developed enough. But for the problem of single 
point source in GR, in its ID formulation, used in this 
article, one can use the classical complex analysis. More- 
over, using the well known complex representation of dis- 
tributions (see Bremermann in it is not difficult to 
generalize the classical results for the case when the dif- 
ferential equation has on the right hand side a distri- 
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bution, which depends on the variable z. Such a term 
may lead to a discontinuity of the solution w{z), or of its 
derivatives. The corresponding complex analysis of the 
solution remains an important open mathematical prob- 
lem, together with the problem of finding and classifying 
all static solutions of Einstein equations with one point 
singularity and spherical symmetry. 

The most important physical problem becomes to find 
criteria for an experimental and/or observational probe 
of different solutions of Einstein equations with spherical 
symmetry. Since there exist physical and geometrical dif- 
ferences between the solutions of this type, one can find 
a real way of distinguishing them experimentally. Es- 
pecially, a problem of the present day is to answer the 
question about which one of the spherically symmetric so- 
lution gives the right description of the observed CDO, if 
one hopes that the CDO can be considered in a good ap- 
proximation as spherically symmetric objects with some 
Keplerian mass M and some mass defect ratio g S (0, 1). 
In this case, according to the Birkhoff theorem, we must 
use the new solutions <|29|) for description of their grav- 
itational field in the outer vacuum domain, outside the 
CDO. 

On the same reason the solutions l|29[l must be used, 
too, for description of the vacuum gravitational field of 
spherically symmetric bodies of finite size, made of usual 
matter, outside their radius R. Because of the specific 
character of the correction of the Newton's low of gravi- 
tation, defined by the term G N Mc~ 2 / In (-^j 1 ) in the po- 
tential (p G (r; M, Mo) 130|) . the difference between our so- 
lutions H29|) and the Hilbert one is in higher order rel- 
ativistic terms, when In ~ 1. This difference will 

not influence the rho-gauge invariant local effects, like 
perihelion shift, deviation of light rays, time-dilation of 
electromagnetic pulses, and so on, but it may be essen- 
tial for the quantities, which depend on the precise form 
of the potential <f G (r; M, Mo) lj3"0")l . For all bodies in the 
solar system the term G N Mc~ 2 / In (-^j-) is very small in 
comparison with their radius R: (p G /2R) Earth ~ 10~ 10 , 
(p G /2R) Q ~ 10~ 6 , etc. Nevertheless, the small dif- 
ferences between the Newton's potential and the mod- 
ified one may be observed in the near future in the ex- 
tremely high precision measurements under programs like 
APOLLO, LATOR, etc. (see for example if we will 
be able to find a proper quantities, which depend on the 
precise form of the potential ip G (r; M, Mo) (|30|1 . A better 
possibility for observation of a deviations from Newton's 
gravitational potential and in four-interval may offer neu- 
tron stars, since for them (p G /2R) nir ~ 0.17 — 0.35. 

The results of the present article may have an impor- 
tant impact on the problem of gravitational collapse in 
GR, too. Up to now most of the known to the author in- 
vestigations presuppose to have as a final state of the col- 
lapsing object the Hilbert solution, or a proper general- 
ization of it with some kind of event horizon, in the spirit 
of the cosmic censorship hypothesis. In the last decade 
more attention was paid to the final states with a necked 
singularities. The existence of a two-parametric class of 



regular static solutions without event horizon opens a 
new perspective for these investigations. Most probably 
the proper understanding of the origin of the mass defect 
of point particle, introduced in the present article, will 
be reached by a correct description of the gravitational 
collapse. 

Without any doubts, the inclusion of the new static 
spherically symmetric regular solutions of Einstein equa- 
tions in the corresponding investigations will open new 
perspectives for a further developments in GR, as well as 
in its modern generalizations. 
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APPENDIX A: DIRECT DERIVATION OF FIELD 
EQUATIONS FROM EINSTEIN EQUATIONS 

Using the normal field variables, introduced in Section 
V, one can write down the nonzero mixed components G^ 
of Einstein tensor in the basic regular gauge (<f(r) = 0) 
as follows: 

p 2cp 1 -4ip 2 , x 

^ = ^=^(2pVi'-2pV' 2 '-(pVi) 2 + (M) 2 + e 2v2 ), 

g; = ((M) 2 - (M) 2 +e 2 ^ 2 ) , 

G B e = G\ = -£!!-^(pV^ + (M) 2 -(M) 2 )- (Al) 

Taking into account that the only nonzero component 
of the energy-momentum tensor is T\ and combining the 
Einstein equations Q , one immediately obtains the field 
equations l|13|) and the constraint (|18fl . A slightly more 



general consideration with an arbitrary rho-gauge fix- 
ing function <p(r) yields the equations and the con- 
straint 

A similar derivation produces the equations (JSJ) , which 
present another version of the field equations in spheri- 
cally symmetrical space-times. In this case: 




APPENDIX B: INVARIANT DELTA FUNCTION 
FOR WEYL SOLUTION 



The function p w (r) = \(^/r/ Pg + \/ Pg/tJ > pa 
which describes the Weyl rho-gauge fixing is an one- 
valued function of r G (0, oo), but its inverse func- 
tion r w (p) is a two- valued one: r w (p) = (p) = 

p(l± y/l — Pc/pj ' P € [pa, 00 )- Obtaining informa- 
tion about distant objects only in optical way, i.e. by 
measuring only the luminosity distance p, one is not 
able to choose between the two branches of this func- 
tion and therefore one has to consider both of them 
as possible results of the observations. Hence, the 3D 
space M.^{g m n(Y)}, which appears in the Weyl gauge, 
has to be considered as a two-sheeted Rlemann sur- 
face with a branch point at the 2D sphere p* = pg 
in M.^y {gmh(P) 0) <t>)} ~ the "observable" space with 
new spherical coordinates p, 9, <fi. Because of the limits 
r+(p) — > oo and r~(p) — > for p — > oo, one can write 
down the invariant <5-function in the case of Weyl gauge 
in the form 

Now, taking into account that the equation 1/ p(r) = 
has two solutions: r = and 1/r = and using the stan- 
dard formula for expansion of the distribution 6(1/ p(r)J , 
we obtain easily the representation l|3(i() . 

APPENDIX C: INDEPENDENT NONZERO 
INVARIANTS OF RIEMANN CURVATURE IN 
NORMAL FIELD COORDINATES 

Using the the representation Q of the metric in nor- 
mal field variables after some algebraic manipulations one 
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obtains the following expressions for the possibly inde- 
pendent nonzero invariants of Riemann curvature tensor 
in the problem at hand: 

R = _ 2e 2(vi - 2v2) [ip'{ - 2E 2 - E 3 ) , 




where E 2 := tp 2 -e 2v2 /p 2 and E 3 := tp^ 2 - p 2 2 + e 2 ^ 2 / p 2 . 

As we see, the invariants w\^ 2 and in our problem 
are real. Hence, in it we may have at most 6 independent 
invariants. 

In addition (idi) 3 = 6 (w 2 ) 2 and the metric JjJJ falls 
into the class II of Petrov classification |24| • As a result 
the number of the independent invariants is at most 5. 

For our purposes it is more suitable to use as an inde- 
pendent invariant 

C := -e 2 <^ 2 (2^+6^1 {<p'x-^)-E 2 -2E^ (C2) 

instead of the invariants w\ and w 2 ■ The invariant C has 
the following advantages: 

i) It is linear with respect of the derivative ip'{. This 
property makes meaningful the expression (|C2|) in the 
cases when tp'{ is a distribution; 

ii) It is a homogeneous function of first degree with 
respect to the Weyl conformal tensor C m „,a^; 

iii) It is proportional to an eigenvalue of a proper ten- 
sor, as pointed, for example, in Landau and Lifshitz 0. 

iv) The difference between the scalar curvature R and 
the invariant ljC2f) yields a new invariant 

D:=R-C=3e 2 '<>^(2v' 1 (ip' 1 -<p 2 ) + £ 2 ) (C3) 

which does not include the derivative <p'[. Due to this 
property the invariant D in the case of regular solutions 
l)22[l . (12(it will not contain a Dirac 5-function. One can 
use the new invariant D as an independent one, instead 
of the invariant C. 

As seen from Eq. I|C1II . the following functions of the 
invariants r\, m 3 and TO5: 

e 2( Pl -2 V2 ) ^// + E ^ = ^(maM-rOA 

e 2( Vl -2 V2 ) ^ _ E2 _ E ^ = ^ (4ri „ TO3 / Wl ) /3 ( C4 ) 

are linear with respect to the derivative <p". Using proper 
linear combinations of these invariants and the scalar 
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curvature R one can easily check that the quantities 
E 2 e 2( Vl -2<p 2 ) ^ and £ 3 e 2( Vl -2 V2 ) are independent invari- 
ants too. 

Thus, our final result is that, in general, for the metric 
iJSJ in the basic regular gauge (p = the Riemann curva- 
ture tensor has the following four independent invariants: 

h := e 2 ^- 2 ^ip'{, 

h ■= e 2 ^ 1 " 2 ^ - e 2 ^/p 2 ) , (C5) 



I 3 := e 2 ^ 1 " 2 ^) {-f 12 + v 12 - e 2 ^ jf) /2, 



These invariants are linear with respect to the second 
derivatives of the functions <fi,2 - a property, which is of 
critical importance when we have to work with distribu- 
tions ip'{ 2 . 
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